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Abstract 

We generalize the construction of the multifractal random walk 
(mrw) due to Bacry, Delour and Muzy to take into account the asym- 
metric character of the financial returns. We show how one can include 
in this class of models the observed correlation between past returns 
and future volatilities, in such a way that the scale invariance prop- 
erties of the mrw are preserved. We compute the leading behaviour 
of ^-moments of the process, that behave as power-laws of the time 
lag with an exponent ( g = p — 2p{p — 1)A 2 for even q = 2p, as in the 
symmetric mrw, and as Cq = P + 1 — 2p 2 X 2 — a (q = 2p + 1), where A 
and a are parameters. We show that this extended model reproduces 
the 'harch' effect or 'causal cascade' reported by some authors. We 
illustrate the usefulness of this 'skewed' MRW by computing the result- 
ing shape of the volatility smiles generated by such a process, that we 
compare to approximate cumulant expansions formulas for the implied 
volatility. A large variety of smile surfaces can be reproduced. 



1 Introduction 

Volatility clustering in financial markets is a well known phenomenon. 
More surprising is the fact that volatility correlations are found to 
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be long-ranged, and cannot be characterized by a single correlation 
time. Recent empirical works Jl], [2[ ||, [|, U suggest that the volatil- 
ity correlation function of various assets (in particular stocks) actually 
decays as a power-law of the time lag, with a rather small exponent. 
Furthermore, the volatility fluctuations are found to be close to (but 
not exactly) log-normal J|, fj|]. These observations have lead several 
authors (8], |9j, || to propose an interesting class of multifractal models, 
where the log-volatility is a Gaussian random variable with a corre- 
lation function that decays in time as a logarithm. An explicit con- 
tinuous time construction of such a process was proposed and studied 
in ||, and was showed to exhibit strict multifractal properties, in the 
sense that the even moments of the log-price difference scale with a 
non trivial power of the time lag (more precise statements will be given 
below). This means in particular that the kurtosis of the process de- 
creases only very slowly with the time lag, in contrast with most simple 
models of stochastic volatility, where the kurtosis drops exponentially 
with time. This model is therefore of interest for option pricing, be- 
cause it is consistent with smiles that flatten only very slowly with 
time [10, 11, 12|. The 'multifractal random walk' (mrw) constructed 
in H is however explicitly symmetric, in the sense that all the odd 
moments of the log-price difference are strictly zero. The same is true 
of the 'cascade' construction of Mandelbrot, Fisher and Calvet || [l3|| 
(see also |14] for an interesting causal extension of this model). On 
the other hand, an important stylized fact of the time series of stocks 
and stock indices is the so called 'leverage effect': past price returns 
and future volatilities are negatively correlated |15|. This effect was 
documented quantitatively in (l5|, |l^, 17 1. This effect is particularly 
strong for stock indices, and induces a significant (negative) skewness 
in the distribution of price returns. The aim of this paper is to study 
a generalization of the mrw which accounts for the leverage effect and 
the corresponding skewness. A particularly interesting class of models 
preserves the multiscaling property of the mrw and extends it to odd 
moments. We discuss several financial applications of these models, in 
particular to option pricing. In the presence of non zero skewness, the 
volatility smile itself becomes skewed. Approximate theories, based on 
a cumulant expansion, are also discussed in this context. 
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2 The Multifractal Random Walk (mrw) 



Let (Xt,t > 0) be a stochastic process, with stationary increments. 
We note by 5gX t the increments Xt+e — X t and by M(q,£) the q — th 
moment E(\5eX t \ g ). (X t ,t > 0) is said to be a fractal process if 

M(q,£) = C q f". (1) 

When the function ( q is linear in q, one speaks of monofractal process. 
This is the case of the brownian motion for which Q q = q/2, and 



more generally of self-similar processes jl8, pTOp , for which the following 
equality in law holds 

buX t = b H 8eX t , 

where b is an arbitrary scaling factor. When the function £ g is non 
linear in q, one speaks of multifractal process. In this empha- 
sized in ||, Eq. (|]) can in fact only hold in a certain 'scaling regime', 
I . The construction of (5|] is based on the following discretized 
process: 



X A t(t) = £ 8XaM 5X At [k] = e At [k]e^ ^ 

k=l 



For financial applications, X At (t) can be thought of as the logarithm of 
the price at time t. The quantities £a*[^] and WAt[fc] are independent 
Gaussian variables with the following covariance structure: 



E(e At \j]e At [k}) = a 2 At5 j:k 
E{uj At \j]ijJ At [k)) = \ 2 \n.p At [\j 



where 



PAt[t] 




when |^| < ^ 
otherwise 



T is the integral time beyond which the multifractal scaling ceases to 
hold. e A t[k] is of mean and for the variance of the process to converge 
when At — > 0, one has to choose: 

E(u At [k\) = -Var(uj At [k\) = -A 2 hi|- 
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Under these conditions, it was shown in |5| that the even moments of 
the process are, in the limit At — > 0, given by: 

/ D \ p-2p(p-l)A 2 

M(2p,£) =K 2 J-\ 

with 

K 2p = (2p- 1)!!((T 2 T) P / du!--- / ^ p TT|^-^r 4A2 . 

./o ./o £J 

The explicit computation of the above integral can be found in [20tl , 
that shows that the moments are finite only for q < q* = 2 + 1/2 A 2 , 
beyond which these moments are infinite for all I, >From the definition 
of Cg, we obtain: 

C g = f(l-(</-2)A 2 ), q = 2p. 

One can also compute the correlation of the local volatility, given by 
a[k] = ae WAt ^. It is easily shown that this correlation decays as a 
power-law of the time lag 

E(a[k] q a[k + £} q ) - E{a[k] q ) 2 oc r** . 
Empirically, A 2 is found to be rather small, in the range 0.02 0.1 

3 A Skewed Mult ifr act al Random Walk 

(SMRW) 

3.1 Definition 

We now generalize the construction of || in order to account for the 
leverage effect, where the volatility is correlated with past price re- 
turns. We first consider the following discretized model (we omit the 
At in e and u for simplicity): 

SX At [k] = e[k] e Q W o[k] = u[k] - K(i, k)e[i], 

i<k 

where K(i, k) is a certain kernel describing how the sign of the return 
at time i affects the (log)-volatility at a later time k. We think of 
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K (i, k) as being positive, and therefore the minus sign accounts for 
the sign of the leverage effect. In order to preserve the multiscaling 
properties of the process, we choose the kernel K(i,j) to decay as a 
power-law: 

Kq 

We now compute exactly the first moments of this process and then 
give an approximate formula (to first order in At) for the higher mo- 
ments. 



3.2 The second moment 

Let t = nAt, with 1 <C n <C N = ^7. We find, using the notation 

(...) = E^(...y. 

n—1 n—1 

2u[i]\ 



(x 2 ) = (C^A*M) a > = E< e W a >< e 

i=0 

n—1 

i=0 
n—1 



i=0 

We now make the following assumptions: 

1 

a > - 
2 

^K^At 1 - 213 <C 1 (2) 

which, in particular, requires that (3 < \ when At — > 0. The conditions 
(|H) ensure that the sum J2 k<i jjzrgfm = Y^=i converges and that 
we can replace, to lowest order in Kq, the exponential terms in the 
expression for (X 2 ), by 1. Under these assumptions, we simply have: 

{X 2 t ) = ah, 

as for the simple Brownian motion, or for the symmetric MRW consid- 
ered in j5|. 
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3.3 The third moment 

The third moment, related to the skewness, explicitly reads: 

0<i,j,k<n 

= 3 e wm 2 ^ +2m ) 

0<i<j<n 

= 3a 2 At E ( e ^ +2 ^)(e[i]e 2K ^^) 

0<i<j<n 

( e T, k<t (K(k,i)+2K(kj))e{k]^ e 2Z l<k<j K(k,j)e[k]^ 

Within the conditions (|2|), the exponential terms (i.e. ( K e K<yk ^ tk )) can 
again be set to 1 to first order, provided the sums converge, which is 

A 2 CT 2 

the case whenever a > 1/2. Using the equality (ee Xe ) = Xa 2 e 2 ; W e 
find 



(X?) ~ 6(a 2 At) 2 E K(i,j)p ^pl 3 

0<i<j<n 

/l m 3A 2 „ o. A 2 ^-^ At At 

- 6 V^A^ T £ (( ,_ i)A|)j(b ._ [ + i)Atp , 



0<i<j<n 

dm 
\u - V \ a+2X2 

0<u<v<t 



-0A^r*At-^* // ** (3) 



One easily proves that, under the condition V2 = a + 2 A 2 < 1, 

dudu t 2 ^ 2 



\ u _ v \a+2\2 (2-v 2 )(\- U2 y 

0<u<v<t 

so that our final result is: 

4 3A 2 

<*?> - - (2 :; )( r:, 2 ) ^ 2 -- ^ 

where \i = a — (3 + 4-. Therefore, the third moment grows as a 
power of the time lag t, with an exponent £3 = 2 — z^- Due to the 
conditions (g), is always positive and therefore (X t 3 ) vanishes in 
the limit At — > 0. In other words, the skewness of the continuous- 
limit process disappears. In practice, however, the elementary time 
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scale At, beyond which price increments are uncorrelated, is equal 
to several seconds, even in futures markets, and to several minutes in 
stocks markets. >Prom a theoretical perspective, however, it would be 
interesting to construct a multifractal process which remains skewed in 
the continuous time limit. We will come back to this point in section 

I 

3.4 Higher order moments 

As q grows, the algebra of the explicit computation of the q — th 
moment rapidly becomes messy, and the exact result cannot be com- 
puted. However, as we show in appendix, one can still compute the 
dominating term in the limit At — * 0. 

3.4.1 Odd moments 

We want to evaluate the generic moment: 

(X? +1 ) = Yl ^ ■ ■ ■ e ^ P +i] ■ ■ ■ ^ 2p+l] )- ( 5 ) 

0<h,...,i2 P +i<n 

The sum over 2p + 1 variables contains many different terms but, as 
we justify in the Appendix, the dominating one is: 

V 2p+ i = £ {^Mhf • • • e [ Jp+1 ] 2 e -W+2-fe]+-+2^ P+ i] ) 
0<ji,...,jp+i<n 

(6) 

where all the variables are paired two by two but one. For small Kq, 
we finally find: 

(A t 2p+1 ) ~ -M 2p+ iK At» t p+l - 2p2>? ~ a , (7) 

where M. 2 p+i is a positive numerical prefactor, and \x was defined 
above. Therefore, we find that all odd moments scale as a power of 
the time lag t, with Q q = p + 1 — 2p 2 X 2 — a (q = 2p + 1). In this sense, 
the model considered is multifractal for odd moments as well. 

3.4.2 Even moments 

We now want to evaluate the generic even moment 

(Xt P )= E (e[h]...e[i2 P ]e*W+-+^) 
0<ii,...,i2p<n 
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Using similar arguments, detailed in the appendix, we can show that 
the dominating term is simply: 

0<ji,-,jp<n 

and that we retrieve, for K small, the result of ||: 

(X 2p ) ~ M 2p t p - 2p(p -V x2 . 

Hence, we find that a power-law leverage correlation function, intro- 
duced as we have proposed, does not affect even moments (to lowest 
order in Kq), but give to the odd moments a multifractal behaviour. 

4 An extension to slowly decaying ker- 
nels 

As we have noticed above, a condition for the convergence of sums 
appearing in the exponentials is a > 1/2. A way to relax this is to add 
a large time cut-off to the power-law kernel K(i,j). More precisely, 
we set: 

K e - r CJ-*) 

where r ~ At/T <C 1. We now consider the case condition a < 1/2, 
but impose: 

T 2a - l Kl&t l - 2 P < 1. (8) 

This ensures that the sum in the exponential terms is convergent and 
remains small. Following the same route as above, the third moment 
reads 



0<u<v<t 

Note that the double integral over u, v converges provided a + 2A 2 < 1. 
Under the condition K<T, the exponential term can be dropped and 
we recover the same result as above, Eq. (f|), with \i = a — j3 + 4-. 
Unfortunately, (||) with V ~ At/T impose that a — j3 > in the limit 




" dudv 

a+2A 2 ' 
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At — ► 0. Therefore, the exponent [i is strictly positive and the skewness 
again disappears in the continuous time limit. On the other hand, the 
possibility of choosing a < 1/2 leads to the interesting possibility of a 
skewness that grows with time. Indeed, the third moment (given by 
Eq. (||)), divided by the third power of the volatility, scales as t 1 / 2 ^ 2 , 
with b*2 = a + 2A 2 . Small values of a therefore leads to a growing 
skewness (at least for t < T), whereas a > 1/2 necessarily leads to a 
skewness that decays with time. 

5 Numerical results 

5.1 Methodology 
5.1.1 Simulation 

In order to generate correlated Gaussian random variables (the u's), 



we follow a quite standard procedure explained in [21|, based on the 
utilization of the Fast Fourier Transform [p2[ . 

The major difficulty is to efficiently compute the skewness convolution 
Yli<k ^(*> We first generate 2n independent Gaussian variables 

Vj[i\, and periodize our series as: 



and we define 

K(i,k) 



V[i] if < % < k 

rj[i + 2n] if k - 2n < i < 



K(i,k) if0<k-i<n 
otherwise 



We can now evaluate the skewness convolution by Y2%=k-2n ^(*> ^) e W' 
again using Fast Fourier Transforms. The use of n 'pseudo' past vari- 
ables and the cut-off imposed to K(i, k) are essential to avoid spurious 
correlations between past volatilities and future returns. The system- 
atic use of the Fast Fourier Transform makes the whole procedure 
quite efficient and allows us to simulate long series of our process (in 
the applications we used series of 2 13 variables but longer series can 
be obtained). 



5.1.2 Evaluation of the moments 

We have evaluated the moments by averaging over a large number of 
realizations of the process. This is not such an easy task, because 
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the quantities we deal with can be small and the convergence can 
be very slow due to the existence of long-range correlations in the 
variables {w[i]}, and the slow power-law decay of the kernel K(i,j). 
In order to eliminate spurious skewness effects and to be more accurate, 
we have systematically averaged over the processes generated by the 
variables {e[i\, uj[i]} and by the variables {— e[i], ui[i]}. In the absence 
of asymmetry (i.e. when K(i,j) = 0), the second process is exactly 
the "mirror image" of the first one and therefore all odd moments are 
strictly zero. Conversely in the presence of asymmetry, the second 
process is no longer the "mirror image" of the first one (see Fig [l]), 
and one finds non zero odd moments. By proceeding this way, we are 
sure that any non zero result can indeed be attributed to the leverage 
kernel. 




Figure 1: A realization of a mrw and its "mirror- image" (left), the same for smrw 
(right). In this case we observe the influence of the kernel K which 
induce asymmetry in the process. 



5.2 Some results 

The results we present are based on 10 000 independent realizations 
of the process. A total of 2 13 variables were used for each series. 
The parameters we used are T = 4000Ai, A = 0.175, corresponding 
to A 2 ~ 0.03. We chose two different values of the kernel exponent: 
a = 0.3 with (3 = a + A 2 /2 and a = 0.6 with (3 = 0.45. The strength 
of the leverage effect is taken to be Kq = 0.1, small enough to trust 
our approximate formulas. We have evaluated the different moments 
for times ranging from t = lOAi to t = 1300Ai, i.e. in a region where 
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At Ct<T. We have plotted the logarithm of the different moments 
with as a function of the logarithm of the time lag. We expect to 
obtain a straight line with a slope given by ( q . We indeed observe a 



-2.5 - 0=0.3 




1 1.5 2 2.5 3 1 1.5 2 2.5 3 

login) login) 

Figure 2: Logarithm (base 10) of the 2 nd (left) and 4 th (right) moments of the 
process, with a = 0.6 (x) and a — 0.3 (A)), as a function of the 
logarithm of the time lag. We have shown for comparison the theoretical 
prediction, which is independent of a. The agreement is extremely good 
over the whole time region, both for the exponent and the prefactor. 



nice powerdaw scaling for the second and fourth moment of the process 
(Fig. ^), with a perfect adequation between the theoretical prediction 
and the simulations. The agreement is again very good for the third 
moment (Fig. |3|), in particular for a = 0.3. There is a systematic 
difference for a = 0.6, that tends to disappear for large time lags. This 
is due to the fact that corrections to our asymptotic formulas tend to 
be stronger when a increases. Indeed, in the above calculations (|H), 
we have set: 



^ At At ff dudv 

o<h<n (C* " (5 " i + 1)^) 2A2 * J J [ 

US«<7<rs 0<u<v<t 



v \a+2\ 2 



Although this is correct in the limit At/t — ► 0, there are corrections of 
the order of (At/t) 1_iy2 , which persist when iaj — ► 1. Let us note that 
this finite time correction is very small for the even moments, since in 
this case V2 is replaced by 4A 2 , which is very small (see Fig. ^|). 
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log(n) 

a.0.6 




1 1.5 2 2.5 3 

log(n) 



Figure 3: Logarithm (base 10) of the 3 rd moment of the processes, for a = 0.3 
(top) and a — 0.6 (bottom), and K = 0.1. The adequation between 
simulations and theory is extremely good for a = 0.3 (top), but less con- 
vincing for a — 0.6 (bottom). This discrepancy can be simply explained 
by the faster convergence of the Riemann sums in the case a — 0.3 than 
in the case a = 0.6 (see Fig. || 



6 Comparison with empirical data 

The comparison between the predictions of the (symmetric) MRW and 
real financial data was discussed in details in (5[] . In particular, the 
way to calibrate the parameters A 2 and T is to study the correlation 
function of the log- volatility, that should behave as the log of the time 
lag up to time T, with a slope which is precisely A 2 . While A 2 is quite 
well determined (its typical value is 0.03), the 'integral time' T is much 
less precisely calibrated (it is found to be « 2 years). The elementary 
time interval At is given by the chosen discretization of the time series, 
provided it is larger than the correlation time of the returns (otherwise 
the assumption that (e[i]e[j]) = Sij is not justified). For example, for 
daily data, At = 1 day. 

As far as skewness is concerned, the following 'leverage' correlation 
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Figure 4: Logarithm (base 10) of the rescaled difference between the measured 3 rd 
moment of the process and the theoretical one ((Xf h ) — (Xf neas ) /(Xf h )) 
as a function of the logarithm of the time lag (n = t/ At). This graph 
confirms that the rescaled difference behaves like a power-law of the 
time lag, with an exponent close to 1 — v%, that depends on a. These 
corrections become large when v% — ► 1 (for example for a — 0.6) and 
explain the discrepancy observed in Fig. || 



function was proposed and studied in |16|]: 

(6X\i]6X\j]*) . 

jL(i,i) = 7J- i<J- 

{{5X[k} 2 )) 2 

Using the above definition we find, in the limit K^a 2 /^ -C 1, the 
following relation: 

that allows in principle to calibrate K(i,j) using empirical data, once 

A 2 and T are fixed. (Note that only the combination KqT~ appears 
in the above formulae, so that the uncertainty on T does not transpire 
on the skewness). The leverage correlation function is however found 
to be, in the case of stock indices, close to a pure exponential in |i — j\, 
with a decay time of ten days or so [jO^], in contrast with Eq. (|9|). 
For individual stocks, however, the decay time is much longer. The 
introduction of well defined time scale would ruin the scaling properties 
of the MRW process. We choose to preserve the scaling properties of the 
process for three reasons: (a) a strict multifractal skewed process is an 
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interesting process in its own right, that may have applications outside 
finance (for example in turbulence J23|]), (b) the scaling properties 
of the SMRW might enable one to find exact analytic formulas for, 
e.g. probability distributions, and (c) the slow decay of the leverage 
function allows one to obtain long time persistent skews, and is useful 
in the context of option pricing (see section || below). In other words, 
even if the historical time series do not exhibit a long ranged leverage 
correlation function, the implied distributions, consistent with option 
smiles, might do so. 



7 Volatility asymmetry 

In p3], an interesting time asymmetry in financial series was detected 



using a wavelet analysis (see also y25| for related work). According to 



these authors, "the volatility at large scales influences causally (in the 
future) the volatility at shorter scales", whereas the converse effect is 
much weaker. This finding is actually deeply related to the leverage 
effect (2^] and our model, with its extended time correlations between 
returns and volatility, should be able to reproduce it, as we show now. 
One can quantify the effect reported in [25, 24] by computing a corre- 



lation between volatilities at different scales n and m, defined as: 

(71 — 1 \ 2 / n+m 

i=0 / \j=n 

n—1 n+m n+m 

i=0 j=n 0<ii<i2<n j=n 



Although this quantity is not exactly the one considered in p4 |, one 
expects that it behaves very similarly. The first term C\ is easy to 
compute and is equal, in the limit At — > 0, to: 



Ci = a 4 T 4 



dudv 

|4A 2 

(r + r j — r — r 



(1-4A 2 )(2-4A 2 )' 

where r = nAt and r' = mAt. This result is symmetric in r and r', 
and therefore cannot explain the asymmetry. The second term, on the 
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other hand, is equal to: 



C 2 = 2a 6 K$T 4x2 At x2+2a - 2 f 3 [ 

J 0<U1<U2<T<V<T+t' 



+ 



(U2 — U\) a (v — Ui) a 

We now change variables and set: 

u\ = 
u 2 = 



(v - U2Y 
du\ du2 dv 



(10) 



(«2 — "Ul) A {V — U\) 2X (v — U2) 



t + yr' 



\2\ 2 



and assume first that t' = er with e <C 1, i.e. we compare future short 
scale volatilities with past large scale volatilities. One finds: 

C 2 ocK$eT 3 - 2a ~ 5x2 . 

Now, the opposite case where r = er' with e <C 1, i.e. where we 
compare future large scale volatilities with past small scale volatilities 
leads to: 

C 2 oc K 2 e 2 - a - x2 r 3-2 Q -5A2_ 

The ratio of the second indicator to the first is therefore equal to 
e i-Q-A ^ which g 0es t zer0 provided that a + A 2 < 1. In this case, the 
asymmetry indeed has the sign found in p^] , that is, the correlation 
between past volatility at short scale and future volatility at large scale 
is indeed weaker than the correlation between past volatility at large 
scale and future volatility at short scale. 



8 An application to option pricing: skewed 
volatility smile 

It is a well known fact that Black-Scholes hypotheses are not satisfied 
in practice. In particular, fat tail effects and volatility clustering are 
absent in the Black-Scholes world. This has important consequences 
on option pricing (particularly for exotic options) and hedging. One 
particularly clear symptom is the so called volatility smile: implied 
volatility from option prices is not constant but varies across both 
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strike and maturity, denning a volatility surface. Although this effect 
is widely known and studied by both academics and practitioners, a 
detailed understanding of the implied volatility surface is still one of 
the major challenge of modern finance. 

Several authors |27, 28, 11, 12] have given simple approximate formulas 



for option prices (or for the resulting volatility smile) when the asset is 
modelled by an arbitrary stochastic process, which capture the basic 
mechanisms that lead to non trivial smiles. The idea is to consider 
that the distribution of the price increments is weakly non Gaussian. 
Using a truncated cumulant expansion [29, 27], one finds corrections to 
the constant volatility case in terms of the skewness and kurtosis of the 



underlying process (see also [30] for equivalent results but in a different 
language). For example, when the log-price difference between times 
and t has a variance at, skewness K^,t and kurtosis K^t, the volatility 
smile vt is approximately given by ||l2||F| 



00 



1 + 



K3,t 

3! 



(2<T t - d t ) 



4! 



1 — d t + 3d t at — 3<7 t 



(11) 



log(S /K) +rt + a 2 t /2 



crt 

where So is the current price of the underlying, ao is (true) volatility, 
K the strike and r the interest rate. In the limit where at, rt <C 1 and 
So — K <C Sq, the above formula boils down to the formula obtained 
in [fill in the context of an additive model: 



dt. 



(To 



1 



«3,t 



K 4,i . 



'-^ck- ^(1-d?) 
3! 4! 



(12) 



Sp-K 
Soa t 



Since our process displays both anomalous skewness and kurtosis, we 
expect to observe this smile effect. The above formula is interesting 
as it gives an explicit expression of the option smiles for the SMRW in 
terms of its known skewness and kurtosis. It is therefore important to 
compare the result ( |Tl| ) with an exact smile obtained by Monte-Carlo 
simulations (FigJEI), in order to test the accuracy of the above formula. 



J The formula given in |T^] is in fact slightly different in that it neglects terms in at 
compared to the moneyness d t , which we have kept below. 
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We model the returns of the underlying by a geometric SMRW: 

S S 

— = fj,St + cr 5X t 



(13) 



where X t is a SMRW. We then evaluate the price of an (european) 
option as an unconditional average of the final payoff (St — K)+ = 
max(S , T — K, 0) over many realizations of our process with fj, = r. This 
assumes that 



the option is hedged using the Black-Scholes A [[1G|, [3jJ, which is 
a good approximation for the minimum variance hedging strate- 
gies. In this case, the 'risk neutral' process is indeed given by 



Eq. (|13|). Other hedging schemes are in fact possible, and would 
lead to small corrections to the option price, but this is beyond 



the aim of the present analysis, which is to test Eq. (11) in the 
framework within which it is established. 

one does not attempt to measure the instantaneous volatility 
and to condition the paths to start from this measured value. 
In other words, we study the unconditional, time independent 
volatility surface. Correspondingly, both the price and the hedge 
are independent of the current level of volatility. 




H T=0.25 (simulatioi 

■: T=0.75 

T=1 

T=0.25 (theory) 
- T=0.5 
T=0.75 



h T=0.25 (simulations) 



Figure 5: Implied volatility vt as a function of the strike K, for different maturities 
(t=0.25, 0.5, 0.75 and 1 year), for a = 0.3 (left) and a = 0.6 (right). We 
have chosen r = 4%, o\ — 20%, A = 0.175 and T ~ 2 years (these values 
are chosen according to |^|); the discretization step is here At — 2~ 8 . 
In both case we clearly see the smile effect. 



In Fig.[| we have drawn the volatility smiles given by the Monte-Carlo 
simulation and the theoretical smile given by formula (|TT1), with the 
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theoretical skewness K^f and kurtosis K^f obtained from our above 
analysis: 

- WooT&tt-* ^ (14) 



(2-u 2 )(l-u 2 ) 



g^4A 2 £-4A 2 

K ^ = (2-4A 2 )(l-4A 2 ) " 3 - (15) 

As already mentioned in section ||, the skewness grows with time when 
a < This can be of great interest for financial applications since 
it is often observed that the skewness of the option smile can persist 
for large maturities. We clearly see in (Fig.||) that the non Gaussian 
nature of our process induces an asymmetric smile and, visually, the 
agreement between simulations and theory is satisfying for maturities 
0.5, 0.75, 1 and strikes between 80 and 120. The case of the shortest 
maturity is less convincing, as expected from a cumulant expansion, 
that is in principle only valid for long enough maturities. In order 
to be more precise we then have performed a parametric fit of these 
curves using Eq. (|il"l) to obtain implied values of the skewness (k™ p ) 
and the kurtosis (k™ p ). We have chosen to fit the curves over the 
whole interval ([75 : 125]) of strikes. It is not obvious that this is 
the best choice because the implied volatility for deep in-the-money 
options (small strike) is known to be very sensitive to errors (because 
of a small Vega). We have also computed the empirical skewness and 
kurtosis from the simulations. We report these results in Tab.[j] and 
Tab]2]. The agreement between the theoretical and implied values is 
only fair, even if the order of magnitude is correct. The implied kurtosis 
is smaller than the theoretical one for small maturities with a better 
agreement for large maturities, whereas the implied skewness tends to 
depart from its theoretical value for large maturities. Note that the 
skewness increases with maturity, at variance with most theoretical 
models in which it decreases as 1/Vt. Note also that for a fixed range 
of strikes, the range of relative moneyness decreases with maturity 
and the skewness becomes more difficult to ascertain precisely. In 
order to track the origin of the observed discrepancy^, we have directly 



2 At least in the regime t ; <C T. For t~^> T , the multifractal features disappear and both 
the skewness and kurtosis vanish, as the process converges to a Gaussian. 

3 Note that the agreement between empirical and theoretical volatility is much better 
(see Fig. |). 
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Table 1: Comparison between implied, theoretical and numerically determined 
cumulants for different maturities, for a = 0.3. The theoretical values 
are computed by formulas (|l4| , |l5| ) . The implied values are obtained by a 
fit of the Monte-Carlo values of volatility using Eq. (11), using all strikes. 
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Table 2: Same as in Tab]l], for a = 0.6 



tested the cumulant expansion on the cumulative distribution of our 
process. We show in Fig. ^ the difference between the numerical 
cumulative distribution of the SMRW and the Gaussian distribution for 
two different time lags, t = 0.25 and t = 1, for a = 0.6. We see that 
the cumulant expansion truncated at second order is only qualitatively 
correct in our case, since we observe systematic deviations between the 
empirical values and the predictions. This is particularly true when 
one uses the exact theoretical cumulants, and suggests that higher 
order cumulants cannot be neglected. This is in fact expected, since 
higher order cumulants, for example even ones, only decay with time 
as K2 P ,t ~ t- 2 P(p- 1 ) A2 for the SMRW, instead of t~ p / 2 for iid returns. For 
example, K6,i/ K 4,t ~ in our case. (Actually, cumulants 
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Figure 6: Difference between the distribution function of the smrw (a = 0.6) 
and the normal distribution, for maturity 0.25 (left) and 1 (right). We 
compare it with the cumulant expansion, evaluated with the implied 
values of the cumulants obtained by fit of the volatility smile and with 
the theoretical values of the cumulants ( |l4] , P"5| ) . The results are similar 
in the case a = 0.3. 



even diverge beyond a certain order p* = 1+1/ 4A 2 , see ) . Of course, 
the fit of the cumulative distribution with the implied values extracted 
from the option smile is better, since these values try to correct the 
inadequacy of the truncated expansion. The conclusion of this study is 
that a simple cumulant expansion for the SMRW only gives a qualitative 
description of option smiles. In order to get quantitative results, one 
should use Monte-Carlo paths. However, the scaling properties of the 
model might allow one to find exact analytic solutions in some cases. 
Finally, we want to show many typical smile shapes can be obtained 
with our process by tuning its degree of asymmetry (measured by the 
parameter Kq). We have shown in (Fig.|^) three smiles obtained for 
different values of Kq: 0.1,0.5, 1. In this example a more pronounced 
skewness can be obtained. The term structure of the skewness and of 
the curvature of the smile can be chosen by playing with the parameters 
a and A 2 . 



9 Conclusion and prospects 

In this paper, we have generalized the construction of the multifractal 
random walk (mrw) of [|| to take into account the asymmetric char- 
acter of the financial returns. Our first motivation was to introduce in 
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Figure 7: Implied volatility v t as a function of the strike, for different values of 
the parameter K , for a — 0.3 (left) and a — 0.6 (right). We have 
chosen r = 4%, <j\ — 20%, A = 0.175. In both case we clearly see the 
growing influence of K$ on the asymmetric shape of the smile. 



this class of models the observed correlation between past returns and 
future volatility, in such a way that the scale invariance properties of 
the mrw are preserved. This correlation leads to a power-law skew- 
ness of the returns, and is in fact deeply connected with the so-called 



"harch" effect, or "causal cascade" observed in [gg, |24j]. Such a process 
can be very useful in finance since it captures most of the stylized facts 
of the returns: non Gaussian shape at short scales, volatility clustering, 
long range dependence, and, with the present work, tunable skewness. 
Possible applications in finance are numerous: pricing of options, risk 
management, etc. As illustrated above, the model considered here 
leads to very versatile shapes of volatility surfaces with anomalous 
term structure. In particular, one can account for smiles with skew- 
ness that is constant or even increases with maturity. This would not 
be possible without a slowly decaying return- volatility (leverage) cor- 
relation. From a more academic point of view, we have not been able 
to construct a skewed multifractal process which remains skewed in 
the continuous time limit. Whether this is a fundamental limitation, 
or that alternative constructions are possible, is an open question. It 
might be possible to marry the retarded volatility model of [16] with 
the mrw and obtain a skewed process in the continuous time limit. 
Another promising path would be to study the precise connections be- 
tween the present model and the recent "causal cascade" construction 
of Calvet and Fisher, that extends the work presented in [§] and might 
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allow to introduce the skew in a satisfactory way. We hope to answer 
these questions in future work [32|. 
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Technical appendix 

The computation of the 2p + 1 th moment of X t involves terms of the 
form 

0<il,...,i q <n 

with 

7* G AT 

i 

Writing more explicitly we can expand these terms as: 

^2 (e[n]^e^ ie [ il ])...(e[g^e^*])(e' n£J [ il ] + '-- + ^ a; ^) J] (e Kie[i] ). 

0<ii,...,i q <n —oo<i<n,i^ii,...,iq 

The exponential terms e[i] 7i e^ e ^ follow from the definition of the ibi 
and Ki are complicated prefactors. The above expression looks quite 
involved but as we search for the leading term, we can use (|2|) in order 
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to simplify it. In particular, to the first order in Kq, (e^ i<Ei ) ~ 1. We 
also have 

(e[f e K > eW ) ~ (e[i]' ri (l + ]C ) e[i]+...) 

(e[i] 7i if 7j is even 

K ) e[i} 1+ ^ if 7i is odd 

Due to (H) , we deduce that the lowest order term in Kq is found when a 
maximum number of 7$ are even, that is all 7« except one (because we 
look at the 2p + 1 th moment) are even. Finally, using basic properties 
of the Gaussian vector, {wj}, we find that: 

\ 7*7*A 2 / j, x A 2 (i£ 7 ?-£ 7 i) 



K ' *L \ (l + \i k -i e \)M 



l<k<i<q 



(l + \i k -ii\)AtJ \AtJ 



We want to know which combination of the {7^} leads to the dominat- 
ing term. One way to do this is to study the scaling with At of the 
term corresponding to a given combination of {74}. Combining the 
previous computations, we find that the exponent of At is: 

At l+(a-/3)+A 2 +p(l+2A 2 )-g-lA 2 £Ll 7? _ 

In the limit At — > 0, only the smallest exponent will contribute. The 
only parameters we can adjust are q and Yll=i l1-> with the constraint 
Yli=i 7i = 2p + 1. The two extrema are: 

• q = 2, 7^ = 2p, 7j 2 = 1 which corresponds to an exponent: 

/i 2 = (a-/3) + y + (p-l)(l-2 P A 2 ). 

• q = p + 1, 7i x = • • ■ = 7j = 2, 7i p+1 = 1 which corresponds to an 
exponent: 

A 2 

//p+i = (a-P) + —. 

The difference ^2 — Mp+i = (p ~~ ~~ 2pA 2 ) is positive as soon as 
p < p* = 2^7. Therefore, we can conclude that the behaviour of the 
first 2p* + 1 moments is dominated by only one term, corresponding 
to pairing all points except one. A similar discussion can be made 
for even moments as well. The above argument is clearer in the MRW 
model 01, where the calculation is easier and can be made without 
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any further assumptions. In this model, the odd moments are zero 
and the even moments have a scaling in At given by the exponent 
p(l + 2X 2 ) — q — 4j- Ya=i if- We immediately see that if q = p and 7$ = 
2, this exponent is 0, which is the condition for the existence of a non 
trivial continuous limit [||. 
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